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20.  ABSTRACT  (cont.) 

tlons  were  made  by  Iterative  approximations  rather  than  using  the  minimized 
total  Gibb's  free  energy  method .  The  difference  in  these  methods  is  examined. 

The  chemistry  of  the  transport  system  has  been  examined  in  stages  in  the 
high  tesiperature  source  region.  Three  models  are  proposed  for  the  flat  profile 
reaction  of  HCf  with  the  arsenic  saturated  source  material  which  is  covered  by 
a  crust  of  GaAs  (or  InP).  In  Model  I  HC1  reacts  only  with  GaAs .  In  Model  II 
a  series  of  reactions  occur  resulting  in  constant  crust  throughout  the  growth 
period.  In  Modal  III  there  is  a  constant  build-up  of  crust  during  the  growth 
period.  Additional  reactions  occur  when  there  is  a  temperature  gradient  in 
the  source.  Actual  experimental  data  for  Ga/As  ratios  corresponds  to  a  combine 
tlon  of  these  models. 

The  effect  of  AsCf^  input,  Hj- inert  gas  carrier  mixtures  and  the  addition 
of  HjO,  0j  and  NH^  were  studied.  The  program  Includes  subroutines  for  inten¬ 
tional  dopants.  For  example,  the  observed  dependency  of  S-doplng  on  AsCl^  input 
can  be  explained  by  the  formation  of  Ga§S.  The  program  also  supplies  Informa¬ 
tion  regarding  supersaturation,  i.e.  nucleatlon,  as  well  as  relative  rates  of 
growth  and  etching  of  the  epl  layer. _ — 
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INTRODUCTION 

Epitaxial  GaAa  ia  used  In  a  variety  of  electronic  devices  (LEDs,  FETs,  IMP ATT* , 
etc.).  The  relatively  vide  band  gap  and  high  electron  nobility  of  GaAa  gives  it 
an  edge  over  silicon  in  microwave  applications  including  solid  state  transmitters 
and  receivers  for  radar  and  communication  systems.  This  requires  materials  with 
less  than  10^  cm’3  impurities  and  mobilities  greater  than  120,000  at  77  K.  High 
purity  epitaxial  GaAs  has  been  prepared  using  the  single  bubbler,  single  boat 

Ga/AsC^/Hj  technique.  Good  results  have  also  been  obtained  by  modifications 
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of  the  basic  technique  using  solid  GaAs  ,  a  two  bubbler  system  ,  flat  furnace 
profiles  ,  H2- inert  gas  mixtures  2,8  and  inert  gases  ®  .  In  each  of  these 
modifications  gallium  is  transported  as  GaCl  together  with  a  mixture  of  Asj 
and  As4  from  a  high  temperature  source  region  via  a  carrier  gas  to  a  low  tempera¬ 
ture  deposition  region  where  epitaxial  GaAs  is  grown  on  a  high  resistance  GaAs 
substrate. 

Unfortunately  the  preparation  of  epitaxial  GaAs  is  still  largely  an  art  and 
consistently  good  material  cannot  be  obtained.  This  results  in  low  yields,  high 
cost  and  components  of  questionable  reliability.  To  overcome  these  problems  the 
Electronic  Research  Branch  of  the  Avionics  Laboratory,  WPAFB  has  initiated  a  pro¬ 
gram  to  re-examine  the  source  and  substrate  materials  and  the  transport  and  growth 
conditions  required  for  the  preparation  of  high  mobility  GaAs  using  a  new  reactor 
system  containing  two  boats  and  three  bubblers  (Fig.  1).  In  this  system  purified 
hydrogen  is  used  as  a  carrier  gas  and  supplied  to  a  manifold  consisting  of  three 
mess  controllers  and  three  arsenic  trichloride  bubblers  maintained  at  constant 
temperature.  The  arsenic  trichloride  from  the  first  bubbler  reacts  with  hydrogen 
gas  in  the  high  temperature  region  of  the  reactor  tube  to  form  hydrogen  chloride 
and  an  equilibrium  mixture  of  HC1  and  AS2  and  Aa^  (equation  1). 

1.  AsC13  +  3/2H2  -  nA«n  3HC1 
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The  hydrogen  chloride  formed  In  reaction  (1)  reacts  with  the  walla  of  the 
silica  glass  reactor  to  fora  chloroallanea  and  with  the  Ga  and/or  GaAs  In  the  boat 
to  fora  GaCl.  Eleaantal  arsenic  also  reacts  with  Ga  In  the  boat  to  fora  a  satu¬ 
rated  phase  Gai.xAsji  **  and  GaAs.  Arsenic  trichloride  fron  the  second  bubbler 
passes  through  a  cracking  furnace  where  reaction  (1)  again  occurs.  The  eleaental 
arsenic  Is  condensed  and  the  purified  HC1  is  allowed  to  pass  over  galllia  in  the 
second  boat  where  reaction  (2)  occurs. 

2.  HC1  +  Ga  GaCl  +  1/21^ 

Finally  arsenic  trichloride  fron  the  third  bubbler  enters  the  high  temperature 
region  of  the  reactor  tube  beyond  the  boat  region.  In  the  low  temperature  region 
epitaxial  layers  of  GaAs  are  formed  on  Cr-,0-  or  Sn-doped  substrates  according  to 
the  overall  reaction 

3.  GaCl  +  l/nAsn  -  1/21^  -  GaAs  +  HC1 

where  Asn  Is  the  equilibrium  mixture  of  As2  and  AS4. 

By  using  the  two  boat,  three  bubbler  system  a  number  of  parameters  can  be 
Investigated:  flow  rates,  Ga/As  and  Ga/Cl  ratloa,  Hj /Inert  gas  ratios,  the 
reactor  teisperature  profile  and  the  source  and  deposition  temperatures.  In  par¬ 
ticular  the  effect  of  these  parameters  on  the  stoichiometry  of  the  epllayers  will 
be  Investigated,  l.e.  what  parameters  control  the  formation  of  n  and  p  type  materials. 
OBJECTIVES 

In  order  to  gain  a  better  understanding  of  the  chemistry  of  the  system  and 
to  contribute  to  the  experimintal  design  of  the  proposed  Investigations,  It  seemed 
worthwhile  to  develop  a  computer  program  which  would  incorporate  all  of  the  fea¬ 
tures  of  the  reactor,  Including  all  of  the  outlined  parameters  and  possible  chemi¬ 
cal  reactions.  The  purpose  of  this  report  is  to  discuss  the  chemistry  of  the 
Ga/AsCl3/H2  and  related  systems,  to  briefly  outline  the  selient  features  of  the 
computer  program  mil  Its  thermodynamic  data  bank  and  to  report  on  soma  of  the  findings 


resulting  from  the  chemical  analysis  and  computer  calculations. 

FINDINGS 

The  Chemistry  of  the  Ga«AsCl-»  System 

In  the  high  temperature  (850°C)  region  of  the  reactor  and  in  the  cracking 
furnace  arsenic  trichloride  and  hydrogen  react  to  form  HC1  and  elemental  arsenic 
(1).  In  addition  the  equilibrium 
A.  As4  ■  2Asj 

is  established.  Next  the  hot  HC1  reacts  with  the  hot  silica  glass  of  the  Inlet 
tube  and  reactor  walls  to  form  a  series  of  ehlorosilanes  (5-8)  S10,  SiCl2  and 
elemental  silicon  vapor,  Si,g. 


5. 

Si02,gl 

+ 

(A-x)HCl 

+  xH2 

-  Si^ClA-x 

+ 

2H20 

6. 

Si02,gl 

+ 

h2 

-  SiO 

+ 

h2o 

7. 

Si02,gl 

+ 

2HC1 

+  «2 

-  SiCl2 

+ 

21^0 

8. 

Si02,gl 

+ 

2HJ 

-  Si.g 

+ 

21^0 

Numerous  reactions  occur  with  the  source  materials.  The  boat  itself  reacts 

12 

with  gallium  (reaction  9)  to  form  volatile  S10  and  Ga20  .  The  ehlorosilanes 
also 

9.  SiOj.gl  +  2  Ga  -  SiO  +  Ga20 

react  with  gallium  forming  additional  silicon  vapor  which  can  dissolve  in  the 

10.  SiHjjCl^^  +  (A-x)Ga  -  Si  +  x/2H2  +  (A-x)GaCl 

gallium  13  '  Elemental  silicon  can  also  react  with  hy  drogen  and  the  HC1  intro¬ 
duced  via  ASCI3  beyoad  the  boat  region  to  reform  ehlorosilanes. 

ehlorosilanes  have  been  uaed  ^  as  sources  of  silicon  in  the  prepratlon  of 
Si-doped  epitaxial  GaAs.  It  is  clear  therefore  that  there  are  ample  sources  of 
unintentional  Si-doping  in  the  Ga/AsCl3/H2  system.  But  it  has  also  been  found 
that  Si-doping  decreases  with  increasing  ASCI3  input  ^  . 

Arsenic  formed  either  by  reaction  (1)  in  an  hydrogen  atmosphere  or  by 
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reaction  (11)  la  an  inert  atmosphere  dissolves  in  gallium  forming  Ga^gAs^ 


11.  AsC13 

+  3Ga 

3GaCl  + 

1/nAsn 

12.  Ga 

+  l/nAan  - 

^l-X^X 

where  X  is  the  atomic  fraction  of  arsenic  in  the  liquids  phase  in  equilibrium  with 
gallium  arsenide.  Unfortunately  this  system  readily  supersaturates,  and  upon 
nucleation  forms  a  crust  of  GaAs^  .  it  is  not  known  as  to  whether  this  crust 
is  completely  uniform  and/or  porous  during  the  growth  period. 

If  the  assumption  that  the  crust  is  uniform  and  nonporous  and  the  temperature 
profile  over  the  source  material  is  flat,  then  the  only  material  exposed  to  the 
gas  phase  is  GaAs  which  reacts  with  HC1  (13). 

13.  GaAs  +  HC1  -  GaCl  +  1/2H2  +  l/nAan 

Reaction  (13)  is  the  reverse  of  reaction  (3). 


If  the  crust  of  GaAs  is  porous,  broken  or  there  is  a  temperature  gradient 
along  the  source  material,  additional  reactions  must  be  considered. 


14.  Ga^xAsy 

+ 

(l-X)HCl 

-  (l-X)GaCl  + 

(1-X)/2H2  +  X/nAsn 

15.  Gai_2x(G*A«)x 

+ 

(1-2X)HC1 

-  (l-2X)GaCl  + 

(1-2X)/2H2  +  X/GaAs 

16.  Ga^_2x (GaAs )x 

+ 

Y/nAsn  + 

(1-2X-Y)HC1  -  (l-2X-Y)CaCl  +  (1-2X-Y)/2H2  +  (X+Y)GaAs 

Reaction  (14)  corresponds  to  a  steady  state  condition  with  respect  to  the  crust, 
i.e.  d(GaAs)/dt  “  0,  and  slightly  Increases  the  amount  of  arsenic  in  the  gas  state. 
Reaction  (15)  corresponds  to  the  build  up  of  crust  material  from  dissolved  arsenic. 
Reaction  (16)  is  a  general  reaction  for  the  simultaneous  formation  of  GaAs  from 
dissolved  arsenic  and  from  arsenic  dissolving  in  the  high  temperature  region  of 
the  boat  and  precipitating  as  GaAs  in  the  low  temperature  region  of  the  boat. 

This  formation  corresponds  to  Shaw's  description  of  crusting  under  a  temperature 

gradient.  Finally  reactions  (14  and  45 )  could  occur  if  gallium  dlffuaes  through 
GaAs. 

Volatile  GaCl  produced  either  by  the  reaction  of  arsenic  trichloride  with 
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gal Hub  or  by  Che  reactions  of  hydrogen  chloride  with  the  Ga-Aa  system  passes  down¬ 
stream  via  the  carrier  gas  to  the  cooler  deposition  region.  Here  the  overall 
reaction  for  the  growth  of  GaAs  is  given  by  equation  (3).  Although  the  isechanlsm 
of  this  reaction  is  not  known,  it  has  been  suggested  that  hydrogen  1?  is  involved 
in  the  desorption  of  chloride  from  absorbed  GaCl  on  the  surface  of  the  growing 
GaAs  in  the  Ga/AsCl3/H2  system  whereas  the  removal  of  chloride  in  the  inert  gas 
systems  J  is  due  to  the  disproportion  reaction. 

17.  3GaCl  -  2Ga, 1  +  GaCl3 

Equally  unclear  is  the  role  of  either  As2  or  As^  on  the  growth  of  the  epi  layer. 

Additional  reactions  included  in  this  study  are  those  related  to  the  inten¬ 
tional  doping  of  GaAs  by  sulfur  using  l^S  as  the  source  of  the  dopant. 


18. 

h2s 

m 

Hj  +  S 

19. 

h2s 

+ 

2GaCl 

2HC1 

+  Ga2S 

20. 

nH2S 

+ 

l/nAsn  - 

nH2 

+  nAaS 

Qualitatively  the  In/PC lj/Hj  and  In/PCl3 /Inert  gas  systems  are  similar  to  the 
corresponding  Ga-As  systems.  Programs  for  these  systems  were  Included  in  the 
computer  program. 

THE  COMPUTER  PROGRAM 

The  computer  program  was  designed  for  use  with  a  Hewlett-Packard  Model  9825T 
Calculator.  The  program  contains  four  main  systems: 

1.  Ga/AsCl3/H2,  inert  gas 

2.  Ga/AsClj /inert  gas 

3.  In/PClj/l^,  inert  gas 

4.  In/PCl3/inert  gas 

Each  system  includes  subroutines  for  determining  the  degree  of  supersaturation  in 
the  deposition  region  and  gives  information  regarding  etching  and  relative  growth 
rates.  Each  system  also  provides  for  the  reactions  involving  silicon  compounds 
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isoelectronic  with  InF  and  their  fomula  weights  differ  by  only  3  mats  units. 

With  the  large  numbers  of  FEF  functions  used,  It  was  convenient  to  reduce  then  to 

a  co anon  polynomial  sod  avoid  Interpolation  during  each  set  up  calculations. 

Errors  in  FEE  and  H°298  are  magnified  in  computing  K^.  The  errors  in  FEF 

are  small,  but  major  diacrepanclea  exist  in  Borne  of  the  enthalpy  values. 

An  important  example  in  the  Ga-Aa  system  is  the  choice  of  the  enthalpy  value 

for  reaction  (4).  Although  more  recent  values  (Table  II  B.)  are  available,  most 

authors  continue  Co  use  the  value  of  H°29g  "  73.5  kcal/mole  calculated  from  the 

20 

data  in  Stull  and  Sinke  .  This  is  the  origin  of  the  common  procedure  of  set¬ 
ting  the  As4  partial  pressure  nearly  equal  to  one  fourth  of  the  input  value  of 
AsCl^  and  the  continued  use  of  the  reaction 

GaAs  +  HC1  -  GaCl  +  I/2H2  +  l/4Aa4 

to  describe  the  transport  and  deposition  reactions  of  gallium  and  arsenic. 

Values  for  the  Aa2/As4  ratio  at  various  temperatures  and  ASCI3  input  are 
listed  in  Table  III.  These  ratios  show  that  at  high  temperatures  A&2  is  equal  to 
or  more  important  than  A s4  whereas  the  latter  is  still  the  predominant  species  at 
lower  temperatures.  Ban  ^  observed  an  As2/As4  ratio  of  1.6  at  850°C  and  a  total 
A8H3  ^nPut  1*6  x  10“2  atm.  This  is  somewhat  higher  than  the  value  of  0.73 
calculated  for  these  conditions.  However  from  his  As2  and  As&  values  H°29g  *46.7 

kcal/mole  which  is  close  to  the  maximum  enthalpy  value  of  48.5  kcal/mole  proposed 

63 

by  Joyce,  Foxon  and  Wood  .  These  authors  also  find  that  only  As2  is  Important 
above  600  K  in  the  MBE  of  GaAs.  In  a  systematic  study  of  the  rate  of  growth  of 
GaAs  Shaw  ^  found  that  below  775-750°C  the  growth  rate  is  kinetically  controlled 
with  an  activation  energy  of  49  kcal/mole.  Thus  it  is  possible  that  at  low  temper¬ 
atures  the  rate  determining  step  in  the  formation  of  GaAs  by  CVE  is  the  conversion 
of  Aa4  to  As2. 

A  second  discrepancy  is  in  the  choice  of  the  enthalpy  of  formation  of  the 
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and  Intentional  dopants.  Variables  include  the  temperatures  of  the  source,  the 
reactor  tube  along  the  furnace  profile  and  the  growth  region,  the  vapor  pressure 
and  flow  rates  at  each  bubbler,  the  composition  of  the  carrier  gas,  the  atomic 
fractions  of  arsenic  and  phosphorus  in  gallium  and  indium,  respectively.  The 
program  also  permits  the  determinations  of  the  concentrations  of  volatile  silicon 
compounds  at  any  point  in  the  reactor  and  the  dopant ‘react ions.  Because  the 
experiments  in  Electronics  Laboratory  are  performed  in  an  open  tube  reactor,  the 
total  pressure  in  the  program  was  arbitrarily  set  at  one  atmosphere  and  the  concen¬ 
trations  of  the  volatile  components  are  calculated  in  atmospheres. 

The  major  assumption  in  this  program  is  that  equilibria  is  established  at 
the  hot  surfaces  in  each  region  of  the  reactor  through  heterogenous  catalysis. 

The  method  of  calculation  involved  the  determination  of  Individual  species  by  the 
method  of  iterative  approximations  rather  than  by  minimizing  the  total  Gibb's 
free  energy  18  .  The  latter  implies  a  closed  system  and  may  lead  to  results  which 
are  inconsistent  with  the  open  tube  experiments. 

In  general  the  calculations  obtained  here  agreed  with  those  reported  In  the 
literature.  Major  differences  reflect  choice  of  thermodynamic  values,  the  reaction 
model  and  the  methods  of  calculation.  These  differences  are  discussed  below. 
Equilibrium  constants  were  calculated  from  the  relationship 

21.  logK  -  <<FEF)  -  H°298/T)/2.303R 

where  R  is  the  gas  constant,  M°29g  is  the  standard  enthalpy  of  formation  at  298  K 
and  FEF  is  the  Gibb's  free  energy  function.  The  latter  is  defined  by 

22.  FEF  -  -  (F  -  H°298)/T 

where  F  is  the  Gibb's  free  energy.  Values  of  FEF  and  H°29g  ware  taken  from  stan¬ 
dard  compilations,  recent  literature  and,  in  some  Instances,  calculated  from 
either  heat  capacity  or  spectroscopic  data  (Table  I).  Because  the  required  data 
for  In?  is  not  available,  the  heat  capacity  data  for  AlSb  was  used.  AlSb  is 
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silicon  dioxide  in  Che  reactor  Cube.  Because  of  Che  persistent  use  of  the  emphanism 
"quartz"  for  silica  glass  most  authors  12,65,66  select  the  incorrect  value  for 
the  enthalpy  of  formation  of  Si(>2  g^.  Table  IV  lists  the  concentrations  of  volatile 
silicon  compounds  formed  in  the  reactions  between  HC1  and  silica  glass  in  the  hot 
reaction  tube  using  different  enthalpy  values.  Typical  source  conditions  (850°C, 

1  atm  H2,  .005^6  atm  AsClj)  were  employed  in  these  calculations.  Because  the 
same  conditions  occur  when  the  AsClj  from  the  growth  and  etch  bubblers  enter  the 
hot  reactor  no  allowance  for  reactions  with  the  source  material  were  made.  Column 
two  and  three  list  the  calculated  concentrations  of  silicon  compounds  based  on 
enthalpy  values  of  -215.6  kcal/mol  (quartz)  and  -202.5  kcal/mol  (silica  glass), 
respectively.  By  using  the  former  enthalpy  value  the  amount  of  the  silicon  com¬ 
pounds  formed  and  the  amount  of  Si-doping  is  underestimated  by  a  factor  of  ten. 

The  effects  of  varying  the  concentrations  of  ASCI3  (and  HCD  and  of  hydrogen 
using  a  f^-inert  carrier  gas  on  the  formation  of  silicon  compounds  has  been 
examined.  In  addition  the  effects  of  H2O,  O2  and  NH3  were  calculated.  Finally 
a  study  was  made  of  the  reactions  of  silicon  compounds  formed  in  the  reactor  with 
Ga,  GaAs  and  Ga^.^s^  Typical  caluclations  showing  the  effect  of  AsCl^  in  1  atm 
H2  at  850°C  are  given  in  Table  V.  Increasing  the  AsClj  input  increases  the  forma¬ 
tion  of  HC1  by  reaction  (1).  This  results  in  an  increase  in  the  concentration  of 
S1C14  and  SiHC^.  The  concentration  of  Sii^C^  is  nearly  constant  while  the  con¬ 
centration  of  SIH3CI  decreases.  The  latter  represents  the  effect  of  the  simultane¬ 
ous  increase  in  H2O  with  an  increase  in  SiCl^.  The  major  effect  of  increasing  con¬ 
centrations  of  HC1  is  the  rapid  decrease  in  the  remaining  silicon  compounds:  SIH^, 
S10  and  elemental  silicon  gas. 

Changing  the  hydrogen  content  of  the  carrier  gas  also  has  a  marked  effect  on 


the  formation  of  silicon  compounds.  Typical  data  are  given  in  Table  VI.  The 
concentrations  of  silicon  compounds  listed  in  colixnn  four  are  the  same  as  those 
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given  in  column  three  of  Table  IV.  The  values  in  column  two  are  for  a  ten  per¬ 
cent  excess  of  hydrogen  over  that  required  for  the  reaction  with  AsCl3 •  As  long 
as  there  is  an  excess  of  hydrogen,  the  concentration  of  SiCl4,  and  therefore  the 
total  amount  of  silicon  compounds  and  water  vapor,  is  very  nearly  constant  over  a 
one  hundred- fold  change  in  the  H2  concentration.  However,  decreasing  the  hydrogen 
concentration  decreases  the  concentration  of  the  remaining  compounds.  Again  the 
effect  is  most  noticeable  with  SiH^,  SiC  and  Si.  In  the  complete  absence  of  hydro¬ 
gen,  i.e.  the  inert  gas  systems,  the  only  volatile  silicon  compound  present  is  SiO 
from  reaction  (9). 

Oxygen  reacts  with  hydrogen  to  form  water.  Thus  either  oxygen  or  water  may 
be  added  along  with  the  carrier  gas  to  repress  the  formation  of  volatile  silicon 
compounds.  There  are  some  practical  considerations  here,  e.g.  avoiding  mixing 
hydrogen  with  oxygen  and  water  with  arsenic  trichloride.  Ammonia  will  react  with 
the  hot  silica  glass  to  form  water  and  silicon  nitride.  Not  only  will  the  water 
23.  3S102  +  ANH3  -  Si3N4  +  6H20 

repress  the  formation  of  the  volatile 

silicon  compounds,  but  the  silicon  nitride  will  coat  the  reactor  and  have  the 
potential  of  rendering  the  reactor  inert  to  further  reactions  with  hydrogen  and 
hydrogen  chloride.  Because  the  sublimation  temperature  of  ammonium  chloride  is 
335°C,  the  reaction  between  ammonia  and  hydrogen  chloride  in  the  hot  reactor  can 
be  ignored. 

An  attempt  was  made  to  develop  models  for  the  reactions  between  HC1  and  the 
source  materials.  Model  I  corresponds  to  reaction  (13)  and  assumes  that  only 
solid  GaAs  reacts  with  HC1.  This  is  the  condition  Shaw  ^  used  in  his  calcula¬ 
tions  for  a  single  phase  system.  Under  these  conditions  the  total  arsenic  in  the 
gas  phase  is  given  by 

As  •  AsC13  +  GaCl 
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where  ASCI3  is  the  input  from  the  bubbler  and  GaCl  is  the  amount  of  gallium 
chloride  formed  in  the  reaction. 

Model  II  corresponds  to  reaction  (14)  and  assumes  a  uniform,  thin  crust  of 
GaAs  throughout  the  entire  growth  period  of  the  epi  layer.  Under  these  steady  state 
conditions,  i.e.  d(GaAs/dt)  =  0,  the  mass  balance  for  arsenic  is  given  by 
As  =  AsCl^  +  xGaCl 

where  x  is  the  atom  fraction  of  dissolved  arsenic  in  the  liquid  phase,  Ga^^s^,. 

Model  III  corresponds  to  reaction  (15)  and  assumes  that  gallium  diffuses 
through  the  crust  and  reacts  with  HC1  and  arsenic  forms  GaAs.  In  this  model  the 
liquid  phase  is  represented  by  Ga^_2x(GaAs)x,  the  crust  thickens  and  the  mass 
balance  for  arsenic  is  given  by 

As  a  AsCl^ 

The  values  for  the  Ga/As  ratios  for  the  three  models  are  listed  in  Table  VII. 

In  making  these  calculations  it  was  assumed  that  equilibria  occurred  between  the 
gas  and  condensed  phases.  A  second  assumption  is  that  the  atom  fraction  of  arsenic 
in  Gai-xAsx  is  constant  and  that  the  liquid  phase  is  in  equilibrium  with  the  crust 
material.  Under  these  conditions  the  free  energy  of  dissolved  GaAs  is  equal  to 
that  of  the  crystalline  material,  the  product  of  the  activities  of  gallium  and 
arsenic  in  the  liquid  phase  is  equal  to  one  and  the  standard  state  of  arsenic  in 
the  liquid  phase  is  the  saturated  solution. 

The  value  of  "0.75"  for  the  Ga/As  ratio  is  due  to  the  fact  that  in  Model  I 
one  atom  of  arsenic  would  enter  the  vapor  phase  for  each  molecule  of  GaCl  formed. 
Because  of  the  small  values  of  the  atom  fraction  of  arsenic  in  gallium  at  850°C  H«67 
the  Ga/As  ratios  in  Models  II  and  III  are  essentially  those  due  to  reaction  (2). 

The  difference  in  the  Ga/As  ratios  for  Models  II  and  III  is  due  to  the  slight 
increase  in  the  amount  of  arsenic  in  the  vapor  phase  in  Model  II  (reaction  14). 

Reaction  (16)  corresponds  to  a  fourth  model  in  which  arsenic  dissolves  in 
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Che  liquid  phase  of  a  two-phase  condensed  system  and  crystallizes  out  as  GaAs. 

The  limiting  condition  for  this  model  is  a  closed  system  in  which  the  arsenic  in 
the  gas  phase  is  in  equilibrium  with  the  condensed  phases.  These  conditions 
qualify  for  use  of  the  minimized  Gibb's  free  energy  technique  used  by  Shaw  ^  in 
making  calculations  for  the  two-phase  system.  They  correspond  to  the  conditions 
of  no  transport  and  no  growth. 

In  addition,  the  silicon  compounds  formed  by  the  reaction  of  HC1  with  the  hot 
reactor  also  react  with  the  source  materials.  The  resulting  concentrations  of 
silicon  compounds  for  different  ASCI3  input  are  listed  in  Table  VIII  in  Table  V1IIA 
Model  I  was  used  as  the  basis  for  the  calculations  while  the  results  using  Model  II 
are  given  in  Table  VIIIB.  Results  using  Model  III  are  not  listed  as  they  are 
essentially  the  same  as  for  Model  II. 

The  results  using  Model  I,  where  the  source  reactant  is  GaAs,  show  that 
with  this  model  the  concentrations  of  the  higher  chlorosilanes  increase  rapidly 
with  increasing  AsCl-j  input  while  the  concentrations  of  SiH^,  SiO  and  Si  decrease 
with  increasing  ASCI3  input.  With  Model  II,  where  no  change  occurs  in  the  crust 
material,  the  concentrations  of  the  chlorosilanes  again  increase  with  increasing 
AsCl^  input.  However  with  this  model  there  is  little  or  no  change  in  the  concentra¬ 
tions  of  Sify,  SiO  and  Si  over  a  twenty  fold  change  in  the  AsCl-j  input.  This  dif¬ 
ference  is  due  to  the  very  small  amount  of  HC1  present  in  Model  II.  In  fact  the 
results  for  the  Si  calculations  using  Model  II  would  suggest  that  increasing  the 
ASCI3  should  have  no  effect  upon  unintentional  Si-doping. 

The  computer  program  is  designed  to  provide  information  regarding  the  rela¬ 
tive  rates  of  growth  and  whether  etching  will  occur.  Growth  of  the  epi  layer 
depends  upon  the  fact  that  the  equilibrium  constant  for  reaction  (3)  increases  with 
decreasing  temperature.  Thus  a  system  saturated  with  respect  to  GaCl  and  Asn  at 
high  temperatures  will  be  supersaturated  at  lower  temperatures  and  will  either 
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form  GaAs  by  nucleaelon  or  deposit  GaAs  on  sn  exposed  GaAs  surface.  For  stoichio¬ 
metric  reasons  the  degree  of  supersaturation  will  be  the  smeller  of  either  the 
GaCl  or  As  concentrations  in  excess  of  that  required  for  equilibrium  at  the  lower 
temperature.  The  calculated  relative  growth  rate  will  also  depend  on  the  model 
used.  For  example,  using  Model  II  etching  should  occur  under  typical  experimental 
conditions  (source  850°C,  substrate  750°C,  1  atm  1^,  .00546  atm  AsCl^  from  the 
growth  bubbler)  when  the  AsClj  from  the  etching  bubbler  reaches  .002  atm.  This  is 
in  good  agreement  with  the  experimental  observations  in  the  Electronics  Laboratory 
and  indicates  that  these  calculations  should  be  useful  in  the  design  of  experiments 
with  the  two-boat,  three  bubbler  system. 

The  computer  program  is  also  designed  to  study  intentional  doping.  In  par¬ 
ticular  the  program  can  be  used  to  study  the  reactions  of  the  dopant  in  the  gas 
phase  and  to  calculate  relative  concentrations  of  reactants.  The  latter  relation¬ 
ships  provide  information  regarding  the  mechanism  of  the  doping  process. 

For  example,  in  one  series  of  experiments  Sirdoped  GaAs  was  grown  on  (100) 

substrates  using  H2S  as  the  source  of  dopant.  Two  sets  of  experiments  were  run: 

the  first  was  under  normal  growth  conditions  while  in  the  second  additional  AsCl^ 

was  Introduced  from  the  etch  bubbler.  The  results  of  these  experiments  are  given 

in  Table  IXA  and  IXB.  In  both  sets  of  experiments  the  ratio  of  nj-ng/^S  was 

constant  over  a  five  fold  change  in  H2S  concentration.  Similar  results  have  been 
68  69 

observed  elsewheres  1  .  However  in  the  second  set  of  experiments  in  the 

presence  of  additional  ASCI3,  the  ratio  of  n^-n^/HjS  has  decreased  by  a  factor  of 

nearly  7.4.  This  effect  of  AsCl^  on  S-doplng  was  also  observed  by  Heyen  and  co- 
69 

workers  who  found  that  the  decrease  was  inversely  proportional  to  the  AsClj 
input.  The  greater  effect  observed  here  is  due  to  the  fact  that  in  this  experi¬ 
ment  the  added  ASCI3  va*  introduced  through  the  etch  inlet  whereas  in  the  experi¬ 
ments  by  Heyen  and  co-workers  the  additional  AsCVj  was  passed  over  the  source 
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material.  In  the  latter  experiments  the  source  material  was  saturated  with  arsenic  at 
880°  and  then  cooled  to  850°  before  beginning  the  growth  experiments  ®  .  Under  these 
conditions  a  considerable  crust  of  GaAs  will  form  corresponding  to  Model  I.  For 
this  reason  Model  I  was  used  in  making  the  calculations  for  HC1  given  in  Table  IX. 

If  S-doplng  is  thermodynamically  controlled,  reactions  (18),  (19),  and  (20) 
must  be  considered.  The  ratios  of  S/H^S,  AsS/HjS  and  Ga2S/H2S  were  calculated  at 
each  H2S  concentration.  In  all  cases  the  ratios  were  constant  over  the  range  of 
H2S  used.  The  S/H2S  and  AsS/^S  ratios  were  independent  of  the  added  AsCl^.  How¬ 
ever,  the  G82S/H2S  ratio  decreased  by  a  factor  of  6.9  which  is  in  good  agreement 
with  the  decrease  in  the  value  of  n^-ng/f^S  in  the  two  sets  of  experiments. 

If  Ga2S  is  responsible  for  sulfur  incorporation  in  the  epi  layer,  then  the 
sulfur  will  be  on  the  arsenic  site,  the  chemisorbed  sulfur  will  be  immediately 
bonded  to  three  gallium  atoms  and  the  sulfur  would  be  expected  to  have  a  high 
sticking  coefficient.  This  is  in  agreement  with  the  high  values  of  S-dopant/H2S 
ratios  observed  in  Table  IX  and  elsewheres  , 

In  the  Ga/HCl/As^/l^  system  it  can  be  easily  shown  that 
Ga2S/GaCl  -  K/C 

where  K  is  a  constant  incorporating  the  equilibrium  constants  for  the  above  reactions, 
the  H2  pressure  and  the  partial  pressure  of  t^S  and  C  is  the  HC1  input.  Thus  in 
this  system  S-doping  should  be  inversely  proportional  to  the  HC1  input.  It  is  pos¬ 
sible  that  the  data  shown  in  Fig.  3  of  reference  (69)  obeys  this  relationship. 
DISCUSSIOH 

When  the  source  material  is  covered  by  a  non-porous,  uniform  crust  of  GaAs, 
only  reaction  (13)  occurs.  When  this  reaction  is  complete,  all  of  the  GaAs  will  be 
removed,  the  liquid  phase,  Ca^.xAsx,  will  be  exposed  and  a  discontinuity  should  occur 
in  the  vapor  phase  ratio  of  Ga/As.  A  discontinuity  has  not  been  observed  either  in 
the  Ga/As  ratio  nor  in  the  growth  rate  nor  character  of  the  epi  layer.  Instead 


either  a  thin  layer  of  GaAs  or  two  phases  are  observed  throughout  the  period  of 
growth  , 

If  the  physical  fact  of  a  thin  crust  is  accepted  as  occurring  throughout  the 
period  of  growth,  then  reaction  (13),  the  reaction  with  GaAs,  and  reaction  (15), 
the  formation  of  GaAs  in  the  boat,  must  occur  simultaneously.  The  overall  reaction 
would  be  given  by  reaction  (14)  corresponding  to  the  steady  state  condition.  A 
requirement  of  the  steady  state  condition  is  that  the  GaAs  crust  on  top  of  the 
liquid  phase  is  sufficiently  porous  that  gallium  can  readily  diffuse  through  the 
solid  layer  and  account  for  97-8  per  cent  of  the  total  gallium  found  in  the  vapor 
phase. 

An  alternate  steady  state  hypothesis  involves  three  steps.  First,  GaAs  reacts 
completely  with  HC1  (or  AsCl^)  exposing  the  liquid  phase.  Second,  the  liquid 
phase  reacts  with  HC1.  Simultaneously  arsenic  in  the  vapor  phase  dissolves  to  form 
a  supersaturated  solution  in  gallium.  In  the  third  step  nucleation  occurs,  a  new 
crust  of  GaAs  forms  and  the  process  is  repeated.  This  would  produce  an  oscillation 
in  the  Ga/As  ratio  in  the  arriving  at  the  surface  of  the  growing  epl  layer.  This 
in  turn  may  effect  the  stoichiometry  of  the  epi  layer  and  its  electrical  properties. 
Shaw  ^  did  observe  an  oscillation  in  the  weight  loss  of  a  small  (100  mg)  drop  of 
gallium  but  he  could  not  observe  a  similar  oscillation  in  the  weight  loss  of  a 
larger  quantity  of  gallium  in  a  boat  nor  did  he  correlate  the  oscillation  in  weight 
loss  of  the  droplet  with  Ga/As  ratios. 

A  third  model  involves  the  two  condensed  phases  produced  by  a  temperature 
gradient  along  the  length  of  the  source  material.  This  model  (reaction  16)  requires 
a  decrease  of  the  arsenic  in  the  vapor  phase  and  a  corresponding  Increase  in  the 
Ga/As  ratio.  This  is  in  agreement  with  the  observation  that  the  cooled  gallium 
will  contain  relatively  massive  crystals  at  one  end  and  a  chin  crust  at  the  other 
and  that  this  observation  can  be  used  diagnostically  to  Indicate  a  temperature  profile. 
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The  program  and  calculations  described  in  the  previous  section  are  based  on 
the  assumption  that  thermodynamic  equilibria  is  obtained.  There  are  few  data  which 
can  be  used  to  test  this  assumption  and,  in  every  case,  non-equilibrium  conditions 
were  reported.  For  example.  Ban  examined  reaction  (2)  using  mass  spectrophotometry 
and  observed  that  the  reaction  was  incomplete  at  useful  flow  rates.  Shaw  ^  found 
that  the  rate  of  growth  of  epi  layers  at  temperatures  above  775°C  decreased  with 
increasing  AsCl^  input.  This  effect  could  be  partly  obviated  by  increasing  the 
area  of  the  source  material.  Nonetheless  above  an  ASCI3  input  of  5xl0“3  atm  growth 
rates  still  decreased  indicating  Inefficient  source  reactions.  DlLorenzo  and  Moore  71 
determined  the  Ga/As  ratio  by  measuring  the  weight  loss  of  source  materials  at  known 
AsCl^  concentrations  under  typical  growth  conditions  and  found  a  ratio  of  1.83-1.92 
over  a  fivefold  range  of  AsCl^  input.  Although  this  value  is  Intermediate  between 
that  of  about  0.75  for  the  reaction  between  HC1  and  GaAs  and  2.83  for  the  reaction 
with  a  uniform  porous  crust  of  constant  thickness,  its  constancy  suggests  that  the 
reaction  is  kinetically  controlled  with  the  slow  step  first  order  with  respect  to 
AsClj . 

The  ASCI3  dependency  of  sulfur  doping  reported  here  and  in  reference 

69  can  be  explained  by  the  formation  of  GajS.  Si-doping  is  also  A8CI3  dependent  13,71 

Recent  experiments  by  Cox  and  DlLorenzo  4  who  used  a  bypass  A8CI3  bubbler,  show  that 

unintentional  Si-doping  depends  only  upon  the  total  AsCl^,  i.e.  chloride,  Introduced 

into  the  system.  These  results  indicate  that  chlorosilanes  do  not  react  with  the 

72 

source  materials.  Instead  they  support  the  observations  regarding  the  dependency 
of  silicon  etching  on  HC1  and  the  hypothesis  that  silicon  is  removed  from  the  sur¬ 
face  of  the  epi  layer  as  S1C12. 

FUTURE  WORK 

A  variety  of  CVE  systems  produce  high  quality  epitaxial  GaAs  and  each  has  its 
own  inherent  problems.  Given  the  widespread  study  and  commercial  application  of 
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the  Ga /AsC^/Hj  technique  it  would  appear  worthwhile  to  continue  an  examination  of 
this  system  and  to  use  it  as  a  model  for  related  systems.  A  large  manber  of  empiri¬ 
cal  observations  have  been  made  with  this  system  on  the  conditions  for  GaAa  growth 
rates.  Yet  few  observations  have  been  made  regarding  the  source  reactions  and  the 
volatile  components  within  the  reactor.  For  example,  Shaw  ^  implies  that  under 
appropriate  conditions  equilibrium  will  be  established.  On  the  other  hand  DlLorenxo 

and  Moore  71  report  that  less  than  two- thirds  of  the  HC1  generated  react  with  the 
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source  material  while  Feng  et  al  state  that  the  amount  of  HC1  leaving  the  boat 
region  is  equal  to  the  amount  of  GaCl  generated. 

The  value  of  the  GaCl/As  ratio  is  a  function  of  the  stoichiometry  and  extent 
of  the  reaction  between  the  HC1  (and  AsClj)  and  the  source  materials.  It  is  neces¬ 
sary  to  standardize  the  boat  conditions  in  order  to  be  able  to  compre  data  obtained 
in  different  laboratories  and  to  determine  directly  the  extent  of  the  source  reac¬ 
tions  and  the  amounts  of  HC1,  GaCl  and  As  reaching  the  substrate.  It  is  only  by 
quantifying  the  conditions  within  the  reactor  that  we  will  be  able  to  determine 
the  rate  laws  for  the  determination  of  growth  and  thereby  a  better  understanding 
of  the  growth  process.  It  is  anticipated  that  the  latter  will  lead  to  the  consistent 
preparation  of  high  purity  epitaxial  GaAs. 

At  high  temperatures  GaCl  7^  and  AS2  7^  can  be  detected  by  vibrational  raman 
spectroscopy  which  has  a  sensitivity  at  the  parts  per  thousand  level  7^  .  It  would 
be  worthwhile  to  initiate  a  series  of  experiments  using  this  analytical  method  to 
determine  the  nature  of  the  source  reaction  and  the  concentrations  of  the  arsenic 
and  gallium  compounds  Involved  in  the  CVD  of  GaAs  using  the  Ga/AsCl^/l^  technique. 
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Table  1 


Enthalpy  Values  and  Sourcea  of  Thermodynamic  Function* 


H°298 

FEP 

Substance 

Value  Ref. 

Source 

Ref. 

As,  ref. 

Tab lea 

20 

A«2 

45. 54 

21 

Tab lea 

20 

A»4 

37. 3A 

21 

22 

AaCl 

28.00 

23 

23 

AsC13 

-64.00 

23 

Spectroscopic  data 

24 

AsS 

25.00 

25 

Spectroscopic  data 

25 

Cl2 

Tables 

26 

Ga,l 

Tables 

20 

Ga,v 

65.00 

19 

Tables 

20 

GaAs 

-21.70 

27 

Heat  capacity 

28 

GaCl 

-20.00 

29 

29 

G&CI3 

-102.00 

29 

29 

Ga20 

-20.7 

30 

30 

13.70 

31 

Spectroscopic  data 

31 

«2 

Tables 

26 

HC1 

•22.06 

26 

Tables 

26 

«20,g 

-57.80 

26 

Tables 

26 

-4.88 

26 

Tables 

26 

in 

Tables 

20 

InCl 

-26.00 

33 

Spectroscopic  data 

32 

InCl3 

-102.25 

33 

34 

InP 

-15.4 

35 

Heat  capacity* 

36 

NR, 

-11.04 

26 

Tables 

26 

°2 

Tables 

26 

P,  red 

Tables 

26 

P2 

42.73 

26 

Tables 

26 

pa 

30.77 

26 

Tables 

26 

pcij 

-61.03 

26 

Tables 

26 

S,rh 

Tables 

26 

25 


H°298 


Substance 

Value 

Ref. 

SI, 8 

Si,g 

107.70 

20 

SiCl4 

-157.10 

26 

S1HC13 

-119.30 

37 

SIHjCIj 

-75.29 

37 

sih3ci 

-32.40 

37 

SiH^ 

8.20 

38 

SiO 

-24.00 

26 

S102,gl 

-202.50 

39 

S13N4 

-127.00 

26 

■Data  for  InP  unaval labia.  Uaad 


FEF 


Source 

Ref. 

Tablea 

26 

Tables 

26 

Tablea 

26 

Tablea 

26 

Tablea 

26 

Tables 

26 

Tables 

26 

Tables 

26 

Tables 

39 

Tables 

26 

heat  capacity  data  for  AlSb. 
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Table  II 

Choice  of  Enthalpy  Values  for  the  Reaction 
As^  *  2A&2 


A*  A H°298 

■  73.5  kc/mole. 

Ref  20 

Year 

Reference 

Reference  Cited 

1964 

40 

20 

1966 

41 

20 

1968 

42 

40 

1970 

43 

40 

1970 

44 

20 

1971 

45 

40 

1971 

10 

44 

1971 

46 

40 

1972 

47 

20 

1974 

9 

44 

1975 

48 

43 

1975 

49 

20 

1975 

50 

43 

1975 

8 

44 

1976 

51 

50 

1976 

52 

48 

1976 

53 

44 

1977 

54 

47 

1978 

55 

44 

1979 

56 

44 

1980 

57 

44 

1980 

58 

56 

1980 

14 

44 

27 


Table  II  (cont.) 

B.  Recent  Experimental  Values 


Year 

Reference 

4H°„ 

JH  kcal/mole 

1967 

59 

62.50 

1973 

60 

54.26 

1974 

61 

54.40 

1975 

21 

53.74 

C.  Additional  Choices  of  Ah°298 

Year 

Reference 

Refei 

rence  Cited 

1975 

17 

59 

1976 

62 

60 

1981 

This  Report 

21 

28 


Table  III 

Asj/As^  Ratios  for  H°29g  ■  53.74 

T°C 


AsC13 

Input, 


Atm. 

700 

210 

800 

850 

900 

.001 

.46 

.973 

2.07 

4.52 

10.15 

.002 

.31 

.65 

1.32 

2.72 

5.72 

.003 

.25 

.52 

1.03 

2.06 

4.18 

.004 

.22 

.44 

.87 

1.70 

3.38 

.005 

.19 

.39 

.76 

1.48 

2.87 

.006 

.175 

.35 

.69 

1.32 

2.53 

.007 

.16 

.33 

.63 

1.20 

2.28 

.008 

.15 

.30 

.58 

1.10 

2.08 

.009 

.14 

.28 

.55 

1.03 

1.92 

.010 

.13 

.27 

.51 

.96 

1.79 

Table  IV 


The  Enthalpy  of  S102  and  Formation  of  Silicon  Com¬ 
pound  a 

Temperature,  850°C 
AsClo  input,  .00546  a  on 

^2*  *tm  Concentration,  Atm 


h°298 


Compound 

-215.6* 

-202. 5b 

SiCl4 

3.0xl0‘7 

2.2xl0“6 

SiHCl3 

7.5xl0-7 

5.3xl0*6 

SiHjClj 

4.0x10“® 

2.9xl0“7 

SiRjCl 

1.2xl0"9 

8.9xl0"9 

SiH4 

5. 2x10“ 11 

3. 7x10* 10 

SiO 

4. 9x10“ 10 

2.5X10"8 

Si 

3. 9x10" 18 

2. 7x10" 17 

h2o 

2.2xl0“8 

1.6xl0“5 

HC1 

1.6xl0~2 

1.6xl0'2 

“Quart* 
bSi02,  glasa 


Tabic  V 


Effect  of  ASCI3  InPuC  on  SI  Compound* 

Temperature,  850° 

H2,  1  atm 

4H°298,  S102  *  -202.5  kcal/mole 


Compound 

0.001 

S1C14 

7.0x10“* 

SIHCI3 

9. 6x10" 7 

S1H2C12 

2. 8x10“ 7 

S1HC1 

4.8x10”* 

Sify 

1.1x10”* 

S10 

1.3xlO“7 

Si 

8.1x10-16 

HjO 

2.9xl0-6 

HC1 

2.99xl0”3 

SI  Compound,  Atm 
ASCI3,  Ata 


0.004 

0.016 

1.2x10-6 

1.4xl0”5 

4.1x10“* 

1.2xl0“5 

3. 0x10“ 7 

2.0xl0"7 

1.3x10”* 

2. 2x10” 9 

7. 0x10” 10 

3.1X10”11 

3.4x10”* 

7.3xl0”9 

5. 3x10“ 17 

2. 4x10” 18 

5.5x10“* 

1.1x10“* 

1.2x10” 2 

4.8xl0”2 

TabU  VI 


The  EffecC  of  H2  on  Si  Compounds 


Temperature, 
AaCl3  Input, 
Ah°298*  Si02 

850°C 

.00546 

-  -202.5  kcai/moie 

Concentration, 

Atm 

Compound 

.00901* 

H2,  Atm 

al 

L& 

sici4 

5.1x10“® 

4.4x10"® 

2.2x10'® 

sihci3 

1.1x10“® 

1.0x10“® 

5.3x10"® 

SiH2Cl2 

4. 6x10" 13 

5. 0x10" 9 

2.9xl0“7 

SiH3Cl 

1.2x10" 17 

1.44X10"11 

8.9xl0"9 

SiH^ 

4.0xl0“22 

5. 5x10" 14 

3. 7xl0‘ 10 

SiO 

3.0X10"11 

3.3xl0-9 

2.5x10"® 

Si 

4. 0x10" 23 

4. 8x10' 19 

2. 7x10" 17 

h2o 

l.OxlO"5 

l.lxlO*5 

1.6xl0"5 

HC1 

1.64xl0"2 

1.64xl0’2 

1.64xl0-2 

*10X  excesa  over  Che  amount  required  to  react  with  ASCI3. 
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Table  VII 


Ga/As  Ratio  as  a  Function  of  AsCl^ 


Temperature,  850°C 
H2,  1  atm 

Ga/As 


AsClj,  atm  Model 

I 

II 

III 

.0005 

0.74 

2.83 

2.99 

.0025 

0.73 

2.83 

2.99 

.01 

0.73 

2.83 

2.99 

.02 

0.73 

2.83 

2.99 

"Completion" 

0.75 

2.91 

3.00 

I 
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Table  VIII 

Effect  of  Source  Material  on  Si  Compound* 

Temperature,  850°C 
H2 ,  1  atm 

dH°29g,  S102  *  -202.5  keal/mole 


A.  Model  I 

AsClj  Input 

Compound 

0.001 

0.005 

0.02 

stci4 

9.2xl0"12 

l.OxlO"8 

4.8xl0"7 

S1HC13 

2.2x10"® 

3  .OxlO"7 

2.6xl0"6 

SiH2Cl2 

l.lxlO"8 

2. OxlO"7 

3. 0x10“ 7 

SIH3CI 

3.4xl0“8 

7.6x10"® 

2. OxlO"8 

SiH^ 

1.3x10" 7 

3.9xl0‘8 

1.8x10"® 

SiO 

4. 7x10" 7 

2.5xlO"7 

5.6xl0"8 

Si 

9. 8x10" 15 

2. 9x10" 15 

1.4x10" 16 

HjO 

8.3xl0"7 

1.5xl0"6 

6.8xl0"6 

HC1 

1.7xl0"4 

1.3xlO"3 

7.4xl0"3 

B.  Model  II 

sici4 

3. 6x10" 17 

6. 1x10" 12 

SlHCl^ 

2. 0x10" 13 

1.6x10"® 

SIHjClj 

2.3xlO-11 

8.9x10"® 

SIH3CI 

1.6x10"® 

2.9xl0"8 

Sify 

1.5xl0"7 

1.3xl0"7 

SiO 

4.9xl0"7 

4.7xl0”7 

Si 

•14 

1.1x10 

l.OxlO"14 

h2o 

7.9xl0"7 

8. OxlO*7 

HC1 

7 .4x10”® 

1.5xl0"4 

Tabic  IX 


Carrir  Cone 

H2S  Cone.,  atm. 
xlO-o 

K.  HC1,  4. 18x10- 3  atm 


.81 

2.0 

1.13 

3.0 

1.58 

3.5 

2.24 

6.0 

3.19 

10.0 

3.8 

8.0 

l.  HC1  -  I.l2xl0’2  atm 


